Draft version February 20, 2009 

Preprint typeset using I4TgX style emulateapj v. 04/20/08 



A SUZAKU SEARCH FOR NON-THERMAL EMISSION AT HARD X-RAY ENERGIES IN THE COMA CLUSTER 
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ABSTRACT 

The brightest cluster radio halo known resides in the Coma cluster of galaxies. The relativistic electrons pro- 
ducing this diffuse synchrotron emission should also produce inverse Compton emission that becomes compet- 
itive with thermal emission fro m the ICM at hard X-ray energies. Thus far, claimed det ections of this emission 
in Coma are controversial (e.g.. lFusco-Femiano et al.ll2004tlRossetti & Molendill2004l) . We present a Suzaku 
HXD-PIN observation of the Coma cluster in order to nail down its non-thermal hard X-ray content. The con- 
tribution of thermal emission to the HXD-PIN spectrum is constrained by simultaneously fitting thermal and 
non- t hermal models to it an d a spatially equivalent spectrum derived from an XMM-Newton mosaic of the Coma 
field dSchuecker et al. 2004). We fail to find statistically significant evidence for non-thermal emission in the 
spectra, which are better described by only a single or multi-temperature model for the ICM. Including system- 
atic uncertainties, we derive a 90% upper limit on the flux of non-thermal emission of 6.0 x 10~ 12 erg s" 1 cm -2 
(20-80 keV, for V = 2.0), which implies a lower limit on the cluster-averaged mag netic field of B > 0.15^tG . 
Our flux upper limit is 2.5 x lower than the detected non-thermal flux from RXTE ( Rep haeli & Gruberil2002l) 
and BeppoSAX dFusco-Femiano et al.l2004l) . However, if the non-thermal hard X-ray emission in Coma is more 
spatially extended then the observed radio halo, the Suzaku HXD-PIN may miss some fraction of the emission. 
A detailed investigati on indicates that ~50-67% of the emission might go unde tected, which could make our 
limit consistent with Rep haeli & Gruberi (120021) and iFusco-Femiano et al.l (l2004|). The thermal in terpretation 
of the hard Coma spectrum is consiste nt with recent analyses of INTEGRAL dEckert et al.ll2007al) and Swift 
(IOkajimaetalJl2 008: Ajel lo et alJl2009h data. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Coma) — intergalactic medium 
— magnetic fields — radiation mechanisms: non-thermal — X-rays: galaxies: clusters 
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1. INTRODUCTION 

In the hierarchical scenario of cosmic structure forma- 
tion, clusters of galaxies form at late times through merg- 
ers between subclusters and through the accretion of galax- 
ies and galaxy groups. The distribution of their massive 
halos in space and time depend sensitively on the under- 
lying cosmology, and much effort has been made to con- 
nect observable properties of the gas to the total cluster 
mass in order to constrain cosmological parameters (e.g., 
Mantz et al. 2008, and references therein). However, merger 
proce sses are known to sign i ficantly disrupt the therm al gas 
(e.g. iRicker & Sarazinll200lt iRitchie & Thomasl 120021) . typ- 
ically biasing inferred masses and the resultant cosmologi - 
cal parameter estimates (Rand aTTet al.l2002tlWik et"aT] 2008). 
Merger-induced shocks and turbulence, besides heating the 
gas, are thought to also re-accelerate relati vistic particles 
present in the intracl uster medium (ICM) dSarazinl 1 19991: 
iBrunetti ife Blasi 2005). Non-thermal electrons, observed via 
diffuse, radio synchrotron emission, have bee n detected in 
over 50 clusters, all of them undergoing mergers (IBuotd2001t 
iSchuecker et al.ll2001l) . If the energy in a relativistic phase of 
the ICM is large enough to add pressure support to the thermal 
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gas, even transiently, the ability to derive masses and there- 
fore use clusters as co smological probes may be compromised 
(Skillman et al. 2008). An assessment of the relativistic con- 
tribution to the energy budget of clusters is necessary to fully 
characterize the state of the ICM. 

Diffuse, cluster- wide synchrotron radio emission, called ra- 
dio halos or relics depending on their morphology, imply that 
both magnetic fields and relativistic electron populations are 
present on large scales. The total luminosity of a synchrotron- 
emitting electron is given by 



(1) 



where oj is the Thomson cross-section, c is the speed of light, 
7 is the Lorentz factor of the electron, and eg = B 2 /%ir is the 
energy density of the magnetic field. For many such electrons, 
the value of Lr depends both on the number of electrons and 
on B and cannot independently determine either. However, 
these same electrons will up-scatter cosmic microwave back- 
ground (CMB) photons through inverse Compton (IC) inter- 
actions, which have a luminosity Lx equivalent in form to 
equation ([]]) but with eb replaced by the energy density of the 
CMB. Since both luminosities are proportional to the number 
of electrons, their ratio gives the volume-averaged magnetic 
field, 

L R _ B 2 /8tt 

L X aT CMB ' 



(2) 



where a is the radiation constant and Tqmb is the tempera- 
ture of the CMB. T he IC radiation should be observable at 
hard X-ray energies (Rephaelii ll977h . Thus far, IC emission 
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has o nly been detected at low significance dNevalainen et al.l 
2004) or, in one case, in a cluster with uncer tain radio emis- 
sion dEckertetal]|2008t but see also lAiello et alj 120091 and 
iFujita et al] 120081) . The measurement of an IC flux from a 
synchrotron source directly leads to a simultaneous determi- 
nation o f the average value of i? and the relativistic electron 
density (iHarris & Romanishinl [19741 ISarazinlfl988l) . There- 
fore searches for IC emission coincident with a radio halo or 
relic are an excellent way to constrain the contribution of rel- 
ativistic materials in clusters. 

The f i rst, an d brightest, radio halo was discovered by 
IWillsonl (119701) in the Coma cluster, and its radio proper- 
ties h a ve perhaps been the best studied (e.g . Giovannini et al. 

2003}. Coma 



[1991 iDeiss et al.l [l997t iThierbach etal. 



has been observed by all the major obse r vatories with 
hard X-ray capabilit i es (lRephaelietal.lll994t iHughesetafl 
Il993t iBazzano et all fl990t Henriksen & Mushotzkyl 119861) . 

and more recently non-therm al detections have been 

claimed by iRephaeli & Gruberl (120 02) with RXTE and by 
iFusco-Femiano et all d 19991 120041) with BeppoSAX, though 
the latter detection is contro versial (Rossetti & Molend il2004t 
iFusco-Femiano et~aDl2007t) . Due to the large field of view 
(FOV) of these non-imaging instruments and the simple char- 
acterization of the thermal gas, the source of this emission re- 
mains uncertain. Even more recently, long (~ 1 Msec) obser- 
vations with INTEGRAL have imaged extended diffuse hard 
X-ray emission from Coma, though it w as found to be com- 
pletely consi stent with thermal emis sion dRenaud et alJl2006c 
lEckert et al]l2007atlLutovinov et al.ll2008l) ~ 

In this study, we present a Suzaku HXD-PIN observation 
of the Coma cluster in an effort to detect non-thermal emis- 
sion associated with the radio halo and potentially confirm the 
RXTE and BeppoSAX detections. The HXD-PIN instrument 
has a non-imaging collimator like those on-board RXTE and 
BeppoSAX, but with a FOV about a quarter as large, which 
reduces the possible contamination from hard point sources 
(Taka hashi et al.l2007l) . Also, the Suzaku particle background 
is ~ 5 x lower than th e backgrounds of either RXTE or Bep- 
poSAX (Mitsuda et al. 2007). In order to clearly distinguish 
the thermal and non-thermal emission visible within the PIN, 
the hard Suzaku PIN spectrum is jointly fit with a spatially 
equivalent XMM-Newton EPIC-pn spectrum. The XMM spec- 
trum, at lower energies and completely dominated by thermal 
emission, allows Coma's thermal and non-thermal properties 
to be simultaneously determined. The XMM and HXD-PIN 
observations are reported in § [2] and the extraction of the re- 
sulting spectra is discussed in § [3] Fits to the joint spectra 
are described in § @] In § [5] we discuss the implications of 
our results for the nature of the hard X-ray emission from the 
Coma cluster. We assume a flat cosmology with J7m = 0.23 
and Hq = 72 km/s/Mpc and a luminosity distance to Coma of 
98.4 Mpc. Unless otherwise stated, all uncertainties are given 
at the 90% confidence level. 

2. OBSERVATIONS 

The Suzaku observation was undertaken as part of AO-1 
from 2006 May 31 through June 4, soon after 16 of the 64 
PIN diode bias voltages were lowered from 500V to 400V, 
but before an additional 16 diodes were similarly lowered. 
We analyze Version 2 of the pre-processed data (PROCVER 
2.0.6.13), which allows for the diode bias drop, with HEA- 
soft 6.4.0 and XSPEC 12.4.0w. For the HXD-PIN instru- 
ment, the standard data selection criteria are applied to ex- 
tract the source spectrum, and the same criteria are used to 



select times for the modeled non-X-ray background (NXB) 
spectrum. Specifically, we select observing times when the 
geomagnetic cut-off rigidity is above a critical value (COR 

> 6 GV), when the satellite is not within the South Atlantic 
Anomaly (SAA_HXD = 0) or has just left it (T_SAA_HXD 

> 500 s), and when Suzaku is pointed above and at least 5 de- 
grees away from the Earth's horizon (ELV > 5°). The strength 
of the NXB is known to be roughly inversely proportional to 
the value of the COR and to be elevated inside the SAA, grad- 
ually decaying to typical levels after SAA passage. These cri- 
teria ensure that the low NXB of the HXD is minimized and 
can be well characterized, which is necessary if it is to be 
accurately modeled. After event selection, the PIN exposure 
time is reduced from 166.2 ks to 156. 1 ks after dead-time cor- 
rection. The HXD-GSO spectrum is found to be consistent 
with the background, so we do not consider it further here. 
We use the Suzaku CCD data from the XIS0 chip to check 
the cross-calibration of Suzaku and XMM-Newton. Standard 
event selection was applied to the XIS0 data, leading to an 
exposure time of 178.7 ks. 

The mosaic XMM-Newton observations of the Coma clus- 
ter, including 14 separate pointings, were done as a part of 
an instrument performance ve rification program, a co mplete 
log of which is presented in Schueck er et alJ ((2004). The 
initial observation s were undertaken by and first reported in 
iBriel et alJ J200l . Seven new observations, aimed at resolv- 
ing the temperature structure of the Coma center, have also 
been performed (PI P. Schuecker). However, high solar ac- 
tivity during the exposures resulted in a high detector back- 
ground above 2 keV, making these observations less suitable 
for our purposes, and therefore we u se only the observation s 
reported in Section 4 and Table 2 of Schuecker et al. (2004). 
We choose only to include the EPIC-pn data from XMM in our 
analysis. Because these observations were made early in the 
mission, they cannot be processed with the standard software, 
though the EPIC-pn data have undergone in-house processing. 
Also, its effective area at high energies is higher than for the 
EPIC-MOS detector, making it the more suitable instrument. 
The benefit of including the EPIC-MOS data is unclear, due to 
the addition of cross-calibration errors and given the already 
high signal-to-noise of the EPIC-pn data. 

3. EXTRACTION AND CONSTRUCTION OF SPECTRA 

To produce complementary spectra from the XMM EPIC-pn 
and Suzaku HXD-PIN data that can be simultaneously fit, the 
background and responses of both instruments must be care- 
fully considered to minimize systematic uncertainties. The 
expected non-thermal signal is near the limit of the PIN sen- 
sitivity, and a robust characterization of this emission partic- 
ularly depends on the accuracy of the PIN background and 
XMM-Suzaku cross-normalization. 

3.1. HXD-PIN Spectrum and Non-X-ray Background 

As the HXD is a non-imaging instrument, we simply extract 
the PIN spectrum from the selected events and group the spec- 
tral bins so that each bin contains at least 30 counts to ensure 
that Gaussian statistics and \ 2 fitting are valid. The response 
matrix is provided in the Suzaku CALDB 8 for Version 2 data 
products, and we use ae_hxd_pinhxnome2_20080129.rsp 
for all source components other than the cos- 
mic X-ray background (CXB), for which 
ae_hxd_pinflate2_20080129.rsp is used. 

8 http://suzaku.gsfc.nasa.gov/docs/heasarc/caldb/suzaku/ 
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Energy (keV) 

FIG. 1 . — Suzaku HXD-PIN NXB model spectrum (red data points) com- 
pared to the Coma cluster data (uncorrected for background). Note that at 
energies above 45 keV, the NXB dominates the data and that deviations of 
the data above the NXB are confined to individual channels that are simply 
statistical fluctuations or are otherwise imperfectly characterized by the NXB 
model. 

The non-X-ray background for a PIN observation is most 
accurately obtained from a model, as opposed to a compara- 
ble blank field observation. This method is motivated by the 
strong dependence of the background count rate and spectral 
shape on the value of the geomagnetic cut-off rigidity (COR) 
and the time since the passage of the satellite through the 
South Atlantic Anomaly (SAA), quantities which vary and 
have a unique distribution for every observation. The model 
matches the distribution of the COR and other parameters of 
the observed data. A model is also required because there 
is no concurrent measurement of the NXB, such as by nod- 
ding between Coma and a blank field. We use the so-called 
"bgd_d" model for Version 2 processed data, which makes 
use of HXD-GSO information as well as the COR and SAA 
values. This NXB model is shown with the PIN data spec- 
trum, uncorrected for the background, in Figure Q] While the 
shape of the NXB is generally well reproduced, the success of 
the model in determining its overall normalization is ±2.3% 
from 15-40 keV and ±4% from 40-70 keV (Mizuno et al., 
Suzaku Memo 2008-03 9 ). These estimates of the systematic 
error in the NXB are extrapolated to the 90% confidence in- 
terval from the la values derived using Earth-occulted data 
in Suzaku Memo 2008-03. We adopt these values (2.3% from 
12-40 keV and 4% from 40-70 keV) as our estimate of the 
90% systematic error in the PIN NXB. To confirm the accu- 
racy of the model background, we extracted events for both 
the data and model for times when the Earth occults the PIN 
FOV (ELV < -5°, all other selection criteria unchanged); the 
resulting spectra are shown in Figure [2] The fractional dif- 
ference between the model NXB count rate and the Earth- 
occulted data is (0.2 ± 1 .0)% over the energies 12-40 keV and 
(-2.0 ±2.6)% from 40-70 keV (lcr errors). Over the whole 
range considered, the fractional difference in count rates is 
extremely small: (0.005 ± 0.9)%. Because of the excellent 
agreement, we do not adjust the level of the background as 
proposed in Ishida et al., Suzaku Memo 2007-10 10 . 

The estimate of the systematic error adopted here comes 
from an analysis of Earth-occulted data, which is the same 
data used to generate models of the NXB. However, it is pos- 

http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2008-03.pdf 

10 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2007-10.pdf 




Energy (keV) 

FIG. 2. — Suzaku HXD-PIN Earth-occulted data (black data points) com- 
pared to the NXB model spectrum for the same time periods of Earfh- 
occultation (red data points). During Earth-occultation, the only events 
should be due to the NXB. Note that the NXB agrees well with the normal- 
ization and shape of the Earth-occulted data at all energies >12 keV. (Only 
this range is used in the spectral fits for Coma.) 

sible that a systematic effect could be present during observa- 
tions of the sky that would not exist during Earth-occulted 
observations, and so it would not be included in the NXB 
model or the estimate of the systematic uncertainty. Mizuno 
et al. (Suzaku Memo 2008-03) attempt to test for this possi- 
bility with "blank sky" observations and find a larger effective 
systematic uncertainty. It is clear that some part of this uncer- 
tainty is due to the fact that the fields aren't entirely "blank" 
and that the source flux will vary field-to-field. Here we re- 
fer to relatively bright sources not considered to be part of the 
more uniform CXB, which has a variance based on the XIS 
sensitivity to point sources that can be taken into account. 
When considering many observations of one region on the 
sky, so that the contamination from sources will vary less, the 
standard deviation drops from 5.8% derived from many fields 
to 5.0%, both of which includes a statistical error of about 
3.3%. While an additional systematic uncertainty, only in ef- 
fect when observing the sky, cannot be ruled out, this drop 
suggests that systematic error estimates derived from sky ob- 
servations are somewhat conservative. Because the contribu- 
tion of contaminating sources to systematic error estimates is 
thus far not well-characterized, we use the value derived from 
Earth-occulted observations throughout to avoid overestimat- 
ing this error. We rely on the assumption that a full accounting 
of contaminating sources would lead to a systematic error es- 
timate similar to our adopted value. However, using the sam- 
ple of 10 ks exposures of all blank sky observations leads to an 
estimate of the NXB systematic error of 4%, after subtracting 
the statistical error an d the expected CXB fluctuation (see e.g. 
iNakazawa et al.ll2008l) . We consider the effect of raisin g the 
systematic error to this higher value (for E < 40 ke V) in § 14.31 
our results and conclusions remain qualitatively unchanged. 

3.2. XMM EPIC-pn Spectrum 

To constrain the thermal contribution to the PIN hard X-ray 
spectrum, it is very helpful to have a spectrum for the same 
region covered by the HXD-PIN FOV but extending to lower 
energies where the thermal emission is completely dominant. 
This low energy spectrum acts as a lever arm on the thermal 
continuum so that the properties of the thermal gas can be 
extracted simultaneously with a potential non-thermal com- 
ponent. Because the ICM of Coma is not isothermal and its 
projected temperature varies across the cluster, a complimen- 
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FIG. 3. — XMM EPIC-pn 2-7.5 keV wavelet smoothed X-ray surface 
brightness image I Schuecker et al. 2004) with contours of constant Suzaku 
HXD-PIN effective area overlaid. The contours are spaced at 10% intervals 
running from 0% to 90% of the effective area for a point source located at the 
instrument center. In between the contours, the regions labeled Rl, R2, etc., 
correspond to those in Table[T]and in the text. 

tary spectrum at softer energies must follow the spatial sensi- 
tivity of Suzaku's HXD. Otherwise, localized regions of even 
slightly hotter gas could mimic the emission of a non-thermal 
source at hard energies. Since the HXD is made up of 64 
individual collimators with optical axes generally aligned to 
within 4' of each other, we approximate the PIN spatial re- 
sponse as a single perfect collimator with a total square FOV 
of D = 65 '5 on a side, 

R con = (D/2 - 9 x )(D/2 - 9 y )/(D/2) 2 , (3) 

where R m \\ is the fraction of the flux detected at angles of 
(9 X , 6y) from the optical axis along the PIN detector axes, 
relative to a point source located at the center of the HXD 
FOV (9 X = y = 0). We have verified that this model fits the 
spatial response of the PIN very well. The complimentary 
XMM spectrum is constructed based on this spatial vignetting 
of the PIN, which is reasonable for our energy range of inter- 
est (< 70 keV). 

In order to build a spectrum that reflects the PIN vignetting 
with good statistics, we extract spectra from 10 regions of 
roughly equal effective area, as shown in Figure [3] The 
boundaries of the regions are spaced at intervals of 10% of 
the PIN sensitivity to a central point source. Because the 
solid angle subtended by a region increases with its distance 
from the cluster center, it turns out that the count rates of each 
of these XMM spectra are comparable. The same response 
matrix is used for all spectra, epn_ef20_sdY7_medium.rsp, 
and the auxiliary response fil es (arfs) for each re gion are gen- 
erated in the standard way dLumb etalJl200l . The back- 
ground spectrum is der ived from the datasets compiled by 
iRead & Pon man (2003), to which we also apply consis- 
tent flare cleaning criteria. Before summing these spectra, 
weighted by the average PIN sensitivity within each region, 
we scale the arfs so they all agree with the central region (RIO) 
arf at 5 keV, while also scaling the exposure times so the flux 
remains unchanged. Similarly normalized arfs are required to 
ensure that the weighted and summed arf will properly repre- 
sent the response of the final summed spectrum. 

We now describe the procedure for creating the summed 
observed rate spectrum, summed background spectrum, and 



corresponding response. Let 0\ be the observed rate spec- 
trum in spectral channel j for region i, and let B'j be the corre- 
sponding background spectrum. We constructed the weighted 
sums 

10 

Oj = J2^0), (4) 
i=i 

10 

i=i 

where w' is the weight of region ; in the PIN spectrum, based 
on the average value of equation (0) inside the region (Ta- 
ble Q}, normalized by the PIN nominal central point source 
sensitivity. Let R l k be the response matrix for region i, de- 
fined such that given a number flux M' k of photons per unit 
area and time at Earth in spectral channel k, R'j k M' k is the num- 
ber of events per second observed in spectral channel j. In 
the nomenclature of X-ray spectral analysis, R'j k is the "rsp" 
file for region i. The corresponding average response matrix, 
Rjk, must be weighted both by the PIN effective area for the 
regions (the weights w,) and by the number of XMM pho- 
tons incident on each area. To include the weighting by the 
incident flux on each area, we formally assume that to first 
order the spectra in all 10 regions are described by models 
M' k that have identical spectral shapes given by nik, but dif- 
fering normalizations N'. That is, M' k = N'nik- In our case, 
we take wij as an APEC model with T = 8.2 keV and abun- 
dance relative to solar of 0.24, which is the best fit model 
to the XMM-Newton EPIC-pn spectrum for the entire Suzaku 
PIN FOV. While the temperatures from the outer 5 regions 
(R5-R1, TableHJ are somewhat inconsistent with our fiducial 
T, the central regions are weighted more strongly and so it is 
more important to accurately match their spectral shape than 
that of the spectra from R5-R1. The average of the tempera- 
tures from R10-R1 is in good agreement with previous mea- 
surements of Coma's global temperature ( Arnau d~et alJl200ll 
IWatanabe et al.lll999h . Also, temperatures in R10, R9, and 
R8 are just consistent at the 90% level with continuum and 
Fe line ratio fits to the XIS dat a extracted from those regions, 
using the method described in Sato et al] d2008l) . The models 
M' k now differ only in overall flux, given by the APEC nor- 
malization N', so each observed spectrum can be described 
as 

O i r B) = Y J R) k M\ = N i Y J R) k mk- (6) 

k k 

Similarly, we define Rjk as 

O j -B j = N^R jk m kl (7) 

k 

where N is the APEC normalization of the summed spectrum. 
Substituting equations (0]l and © into equation © yields, af- 
ter some algebra, 

k k i 

so it is clear that 

10 N' 

%=E w V^- (9) 

i=l 
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TABLE 1 
XMM Regions and Spectral Fits 







kT 


Norm. a 


Region 


PIN Weight 


(keV) 


(10~ 2 cnT 5 ) 


RIO 


0.933 


8.25±0.21 


1.67 ±0.01 


R9 


0.844 


8.33±0.14 


4.25 ± 0.02 


R8 


0.746 


8.27±0.15 


5.17 ±0.03 


R7 


0.647 


8.07±0.17 


5.09 ±0.03 


R6 


0.547 


8.07 ±0.23 


4.52 ±0.04 


R5 


0.448 


7.40 ±0.34 


3.65 ±0.05 


R4 


0.348 


7.39 ±0.46 


2.93 ±0.05 


R3 


0.248 


6.99 ±0.56 


2.84 ±0.06 


R2 


0.147 


7.65 ± 0.44 


3.30 ±0.05 


Rl 


0.0421 


7.45 ±0.68 


2.76 ±0.06 



11 Normalization of the APEC thermal spectrum, which 
is given by {10~ 14 /[47r(l +z) 2 D\]} J n e n H dV, where 
Z is the redshift, is the angular diameter distance, 
n e is the electron density, tin is the ionized hydrogen 
density, and V is the volume of the cluster. 

The value of the weighted normalization N is given by N = 

For all fits of the XMM spectra, the energy range considered 
differed slightly from the nominal 2-12 keV energy range due 
to calibration issues. At energies near ~ 2 keV, there exists 
a sharp edge in the response due to gold in the mirrors, and 
between 8 and 9 keV there are variable background lines due 
to copper and zinc. Neither of these features can be suffi- 
ciently accounted for given the current calibration, and they 
tend to become important in regions with very good statistics 
(all of the Coma XMM spectra) and/or regions of low surface 
brightness particularly near the outer edges of the detector. 
We avoid these problems by excluding these features and only 
fitting over the range 2.3 < E < 7.5 keV and 9.5 < E < 12 
keV. Also, the gain can vary by up to 30 eV, especially during 
the period when many of these observations took place (Mar- 
cus Kirsch, XMM Calibration Document 11 ) and the redshift 
found from fits differs significantly from the nominal value of 
Z = 0.023. We fit the gain with a linear function, assuming the 
redshift to be 0.0232, using the GAIN command in XSPEC 
for all 10 spectra. We adjusted the gain such that the new 
response energies E' are related to the original energies by 
E' = £71.00519 + 0.010312. None of these calibration issues 
have any important effect on the characterization of the con- 
tinuum features (IC and hard thermal bremsstrahlung) which 
are the subject of this paper. However, not including these ef- 
fects would result in high values of \ 2 f° r the fits due to the 
very good statistics in the Coma XMM EPIC-pn spectra, and 
thus make it more difficult to determine the uncertainties in 
parameters. 

Even after these adjustments, fits to the weighted and 
summed final spectrum with any model produce high chi- 
squared values (x 2 ed > 1.3). A close inspection of the contin- 
uum at various energies reveals that the residuals are slightly 
larger relative to the errors than would be expected by % 2 
statistics, indicating that we have underestimated the errors. 
Because of the high signal-to-noise of the Coma observa- 
tions, the statistical errors no longer completely dominate 
over channel-to-channel systematic effects, caused by, e.g., 
differing/variable charge transfer inefficiencies across the de- 
tectors and/or between observations, which were obtained 
over a two year period. We add a 3% systematic error to the 

1 1 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0018.pdf 



count rate of each channel in the final spectrum to account for 
these uncertainties, which leads to more reasonable values of 

3.3. XMM EPIC-pn and Suzaku HXD-PIN Cross-Calibration 

Joint fits between data from different instruments/missions 
require a careful consideration of the relative overall calibra- 
tion if the validity of fits are to be believed. Instead of di- 
rectly finding the cross-normalization through other obser- 
vations of a spectrally simple source, such as the Crab, we 
choose to use the XIS0 chip of the Suzaku XRT as an interme- 
diary. An advantage of this method is that it does not require 
any assumptions about the stability of the absolute calibra- 
tion of each instrument between calibration observations and 
our observations. To justify our use of the XIS0 data to cal- 
ibrate the absolute flux level, we compare the flux observed 
by XIS0 to the ROSAT 0.5-2 keV flux, which was derived 
using Snowden's ESAS software package dSnowden & Kuntzl 
1998); these fluxes agree to within 1%. 

The overlapping spatial and spectral coverage of the XMM 
EPIC-pn and XIS0 instruments allows a trivial comparison 
of the flux for a region on the sky. We extract an XIS0 
spectrum from the same region as XMM spectrum R10, and 
we generate rmf and arf files for this region using the 2-7.5 
keV wa velet-smoothed image created from the XMM EPIC- 
pn data dSchuecker et al .1120041) . Though the large XIS PSF 
(~ 2') will scatter photons both into and out of this region to 
a much greater degree than occurs for XMM, this effect is ac- 
counted for in the arf and tied to the XMM data. So while 
spatial inhomogeneities will not impact the comparison, the 
shapes of the spectra will not necessarily be identical. The 
overall flux, however, is not sensitive to small variations in 
the temperature, and so it provides a good quantity to estab- 
lish the XMM-XIS cross-normalization. For this region, we 
find that the XMM flux is 15% below the XIS0 flux from 2- 
7.5 keV, and the XMM flux, extrapolated to 0.5-2 keV, is simi- 
larly 15% below the ROSAT flux, and so we scale the summed 
XMM arf by this factor. 

The cross-normalization between the XIS chips and the 
HXD-PIN has been well studied for observations of the Crab 
nebula (Ishida et al., Suzaku Memo 2007-1 1 12 ). We adopt 
their PIN/XIS0 relative normalization factor of 1 . 1 32 ± 0.014, 
increasing the PIN arf, and thus lowering the measured flux, 
by 13.2%. 

The associated systematic error for both cross- 
normalization corrections is estimated to be 1-2%. However, 
the normalization of the R10 spectrum may differ from that 
of the other XMM region spectra, and also the XIS0-PIN 
relative normalization may be different due to the fact that 
Coma is spatially diffuse while the Crab nebula is comparable 
in size to the XIS resolution. These issues suggest that the 
true cross-normalization systematic uncertainty is probably 
larger. We therefore take the combined cross-normalization 
systematic error to be 5%, which is about as large as can be 
reasonably allowed by the simple constraint that a model 
can be continuously fit across the 12 keV boundary between 
the XMM and PIN spectra. Specifically, we vary the cross- 
normalization until the average of the highest signal-to-noise 
PIN channels, covering 12 keV < E < 16 keV, disagrees with 
the model by ~ 2 - 3a. 

12 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2008-03.pdf 
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TABLE 2 
Fits to Joint XMM-PIN spectra 





kT 


Norm." 


T or kT 






Model 


(keV) 


(cm" 5 ) 


( b ) 


Norm. c 


X 2 /dof 


Single T 


8.45 ±0.06 


0.218±0.001 






1676.05/1689 


T+IC d 


8.42 ±0.06 


0.218 ±0.001 


-1.6 


(4.6 ±3.5) X Kr 9 


1671.29/1688 


T+IC d 


8.45 ±0.07 


0.217 ±0.002 


2.0 


(2 + ' 2 ) x 10- 4 


1676.18/1688 


2T e 


8.0 


0.17 


10.1 


0.05 


1672.34/1687 












1684.35/1690 



a See the note following TablefT] 

b Value is Y for the T+IC model and kT (in keV) for the 2T model. 

c Value is the normalization of the power-law component for the T+IC model, which is 

the photon flux at a photon energy of 1 keV in units of photons cm -2 s _1 keV -1 . For the 

2T model, the value is the normalization of the second APEC thermal model (see the note 

following Table [T} in units of cm -5 . 

d Value of r is fixed when deriving errors. 

e Parameters unconstrained. 



3.4. Cosmic X-ray Background 

We modeled t he cosmic X- ray background (CXB) spectrum 
shape following Bol dtld 19871). specifica lly using the analytical 
form proposed by Grube r et al.1 (Q999) based on the HEAO- 
7/A2+A4 data. This shape is well-established over the energy 
range 3 < E < 60 ke V and has been confirm e d in subsequent 
measurements (e.g. IRevnivtsev et all 120031: iChurazov et al.l 
I2007t lAiello et al.ll2008l) . We adopt a 10% larger normaliza- 
tion of the spectrum, relative to the original HEAO-1 determi- 
nation, to agree with the more recent measurements by INTE- 
GRAL (Churazov et al.||2007). This increase is further justi- 
fied by, and consisten t with, the (8 ± 3) % higher normaliza- 
tion found with Swift (lAjello et al.l l2008). Though these most 
recent measurements lie slightly, but systemat i cally , above the 
canonical spectrum, as noted by lAiello et"ail (12008b they are 
not inconsistent with other observations at E > 10 keV. At 
the peak of the CX B spectrum, the me asurement precision 
of HEAO-1 is 10% (Marshall et al.ll 19801) . and the measure- 
ment made with the BeppoSAX PDS is consistent at the 90% 
level with a normalization 12% larger (Fro ntera et al.ll2007h . 
In XSPEC, we model the CXB as 

/ F \" L29 

CXB(£)= 1.056 x 10" 2 — e -£/(4i.i3keV) 
\ 1 keV / 

photons cm" 2 s" 1 keV" 1 , (10) 

where the normalization is set by a 2° x 2° s olid a ngle of the 
sky to match the provided response file (see § 13. It . 

Cosmic variance due to large scale structure depends on the 
solid angle of the observation (57 = 0.32 deg 2 for the PIN) and 
on the cut-off flux of removed point sour ces (5 cut ), determin ed 
by the XMM source completeness (Finoguenov et al. 2004) to 
be S cut (12-70keV) = 2.2 x 10" 13 erg s" 1 cm" 2 . Since the vari- 
ance ctcxb/Icxb fr ' 5 S° u 25 , we can estimate the variance in 
our observation relative to another measurement assuming a 
log/y-lo gS relation of N(S) oc S" 15 . Using the HEAO-1 A2 
estim ate (|Shaferll 19831: iBarcons et alll2000t IRevnivtsev et all 
12001 with n = 15.8deg 2 , S cut = 8 x 10" n erg s" 1 cm" 2 , and 
o~cxb/Icxb = 2.8% (la), we find a variance of 7.4% (90% con- 
fidence), which we take as an additional systematic error in 
the PIN flux. To account for th e 10% discrepancy b etween the 
HEAO-l and the IN TEGRAL (IChurazov et al.l2007l) and Swift 
( lAiello et al.l 120081) observations, we also estimate the stan- 
dard deviation of these measurements to be 7% (90% confi- 



dence). Adding these uncertainties in quadrature, we take the 
total systematic error in the CXB normalization to be 10%. 
Below 20 keV, the CXB emission is <; 10% of the total flux, 
and it just becomes the dominant source of emission at ~ 50 
keV. 

For the XMM data, the background spectra include unre- 
solved point sources that make up the CXB, so they do not 
need to be modeled. 

3.5. Point Sources 

Point sources in the X MM-Newton Com a mosa ic have al- 
ready been identified by Fino guenov et al.l ((2004), who also 
give their count rates in three energy bands (0.5-1 keV, 1-2 
keV, and 2-4.5 keV). For each of the 72 sources, we as- 
sume the spectrum to be described by an unabsorbed power 
law and fit this model to each spectrum. We found that the 
sum of all these models, weighted by w' according to their 
positions, could be more concisely described by the sum of 
two power laws with photon indices 2.1 and 1.6 and normal- 
izations 8.54 x 10" 5 and 1.23 x 10" 4 photons cm" 2 s" 1 keV" 1 
at 1 keV, respectively. While a simple power law description 
poorly characterizes some of the sources, care is taken to en- 
sure that individual fits, when extrapolated to high energies, 
are not unphysical. Their composite spectrum accounts for 
<; 1 % of the XMM flux and is therefore unimportant relative 
to other systematic effects. For this reason we do not go to the 
extra effort to exclude the sources from the XMM spectra. As- 
suming the spectral fits are reasonably valid, the point sources 
account for ~ 10% of the CXB at PIN energies. During fits of 
the joint XMM and Suzaku data, we include this point source 
composite model for both spectra. 

The brightest of these point sources is X Comae, a back- 
ground AGN with a flux ~ 10 x brighter than any other source 
in the field. From XMM RGS observations, it is known to have 
a steep spectrum (T ~ 2.4) and t o vary in flux by a bout a fac- 
tor of 2 over the course of 1 year dTakei et al.l20 07). However, 
due to its position, nearly 90% of the flux from X Comae is 
not detected by the HXD, so this source does not significantly 
contribute to the PIN spectrum. 

4. SPECTRAL FITS 

In our spectral fits, all model components are absorbed by 
the neutral hydrogen column density toward Coma, Nh = 9 x 
10 19 cm 2 [a verage of value s deriv ed from lDickev & Lockmanl 
(1990) and Kal berla et ail d2005l) l. though this absorption is 



Non-Thermal Emission in the Coma Cluster 



7 



negligible at energies above 2 keV. In general, we character- 
ize the dominant thermal emission in XSPEC with the APEC 
model for E < 40 keV and with the MeKa model for E > 40 
keV. In the currect version of XSPEC, the APEC and MeKaL 
models are undefined above 50 keV, though the MeKa and 
bremsstrahlung models are defined. We tie the parameters 
of the MeKa model to the APEC parameters, except for the 
MeKa normalization, which we reduce relative to the APEC 
normalization by 5.5% to bring the models into agreement 
at high energies. Also, the abundances relative to solar and 
t he re dshift are fixed, to 0.24 and 0.0232 respectively (see 
§ |3.21 i. This value for the abundance is based on fits to the final 
weighted and summed XMM spectrum alone, and the best-fit 
abundances of all the individual spectra from the 10 regions 
is also consistent with this value. Since we are interested in 
continuum features, the exact choice for the abundance does 
not strongly affect the results. The spectral fitting results are 
summarized in Table [2] 

4.1. Joint XMM-Newton and Suzaku Spectral Fits Without 
Considering Systematic Errors 

We simultaneously fit the Suzaku HXD-PIN and XMM- 
Newton EPIC-pn spectra for the PIN FOV. First, we consider 
only a single temperature fit, in order to establish whether the 
addition of a non-thermal component actually improves the fit 
(Fig. gj. We find a best fit temperature of 8.45 ± 0.06 keV, 
which is in general agreement with similar fits to the PIN data 
(8.3 ±0.3 keV) and XMM data (8.37 ±0.12 keV) individu- 
ally. Note that the dip at 15 keV is a known problem with 
the current NXB model (Mizuno et al., Suzaku Memo 2008- 
03 13 ). Since each spectrum is individually described by the 
same average temperature, the existence of excess emission 
at hard energies is unlikely. While all of these temperatures 
are slightly higher than the clu ster- wide average temperature 
of 8.2 keV (Hu ghes et al.ll 19931) . the energy range in this and 
similar fits typically extends to energies below 2 keV and thus 
includes more low temperature gas. 

The addition of a power-law non-thermal component pro- 
duces a formally better description of the spectra, according 
to the f-test, improving the overall fit (Table |2), but only for a 
photon index T < 0. Allowing the temperature and power law 
photon index to vary along with each component's normaliza- 
tion, we find T = 8.42 ± 0.06 keV and T = -1.6, though T is 
poorly constrained. If we fix T to this best-fit value, the IC 
component is significant at the 2.2cr level without considering 
the effect of systematic uncertainties. However, this photon 
index is completely i nconsistent with the sp ectral index of the 
radio halo (r ;> 1 .8. lGiovannini et al.l["l993l) . While we might 
expect a flatter spectrum for IC emission, since the hard X-ray 
photons are emitted by somewhat lower energy electrons than 
the radio e mission, and the radio spectrum flattens at lower 
frequencies (Thierb ach et al.l2003l) . a rising IC spectrum with 
energy is completely unexpected and unphysical. The power 
law fit, in contrast to finding an actual power law signature 
in the data, is instead compensating for a slight excess at high 
energies while minimizing its impact on the overall fit at lower 
energies. Notice that the residuals in Figure [4] above 40 keV 
lie systematically, if not significantly, above the model. This 
excess at energies above 40 keV can be explained as a ~ 2% 
underestimate o f the NXB, as suggested by the Earth-occulted 
spectrum (see § 13. It . Increasing the background level by 2% 

13 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2008-03.pdf 
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FIG. 4. — Suzaku HXD-PIN spectrum (E > 12 keV) and the combined 
XMM spectrum (E < 12 keV) corresponding to the spatial sensitivity of the 
PIN. Shown as solid lines are the best fit models for a single temperature 
thermal component. The thermal model ("APEC", green) is nearly coincident 
with the data, though falling below it at higher energies. Also included for all 
joint fits are the the total spectrum for the "XMM Point Sources" (red) and the 
Cosmic X-ray Background ("CXB", purple), the latter of which only applies 
to the PIN spectrum since the CXB is subtracted from the XMM data along 
with the NXB. 
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FIG. 5. — Suzaku HXD-PIN spectrum (E > 12 keV) and the combined 
XMM spectrum (E < 12 keV) corresponding to the spatial sensitivity of the 
PIN. Shown as solid lines are the best fit models for a single temperature ther- 
mal component plus a non-thermal component. The thermal model ("APEC", 
green) is nearly coincident with the data, though falling below it at higher 
energies. The non-thermal model ("Power Law", light blue) is the faintest 
model component for both spectra, and the photon index is fixed at T = 2.0. 
The other two components are described in Figure|4] 

for E > 40 keV results in a best-fit power law component very 
similar to the model used for the XMM point sources, with 
r ~ 1.6, but it is not significant at even the 1<t level. A sim- 
ilar result is found if the photon index is fixed at T = 2 and 
the NXB above 40 keV is not increased; this fit is shown in 
Figure [5] In this case, the fit is not improved by the addition 
of a power-law component to the model. 

Interestingly, a two-temperature model for the ICM yields 
only a slightly better fit to the data than does the single tem- 
perature model (see Fig. [6j, though the addition of a sec- 
ond temperature component is probably not formally justi- 
fied. This result is mainly due to the addition of the 3% sys- 
tematic error to the XMM spectrum. Without including that 
error, a two-temperature model produces a clearly improved 
fit over a single temperature model, indicating that the addi- 
tion of this error is somewhat obscuring evidence for a multi- 
temperature continuum. In either case, the two temperatures 
are not strongly constrained, but they are broadly consis tent 
with the spatial variations in Coma's temperature (see § 14.21 
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FIG. 6. — Suzaku HXD-PIN spectrum (E > 12 keV) and the combined 
XMM spectrum (E < 12 keV) corresponding to the spatial sensitivity of the 
PIN. Shown as solid lines are the best fit models for a two-temperature ther- 
mal component. The thermal model ("APEC+APEC", green) is nearly coin- 
cident with the data, though falling below it at higher energies. The other two 
components are described in Figure[4] 

and Fig. [7J». Therefore, even before systematic errors are con- 
sidered, the case for the inclusion of a non-thermal component 
is not strongly motivated. 
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FIG. 7. — XMM-Newton temperature map across Coma with HXD-PIN con- 
tours of constant PIN effective area overlaid at 10% intervals. The XMM- 
Newton spectra were fit in square spatial regions 4' 3 on a side. The temper- 
atures, given in keV by the color bar, are accurate to either a few tenths of a 
keV (in the center) or 1-2 keV in lower surface brightness regions. Temper- 
atures shown here were determined from fits to the 0.5-14 keV spectrum in 
each region. 



4.2. Multiple Thermal Components in Coma 

In most previous attempts to measure a non-thermal com- 
ponent in the hard spectrum of Coma, the thermal emis- 
sion was modeled as a single temperature plasma character- 
istic of the average global state of the ICM. Ho wever, Coma 
is kn own to host temperature variations (e.g. iHonda et al.l 
119981). Generally, all c lusters exhibit a multi-temperature ICM 
dCavagnolo et alj 120081) . and this is especially true of merg- 
ing cluster s like Coma, which ten d to host hot regions due 
to shocks ( Mar kevitch et al.l [T998). At hard energies, where 
the exponential turnover in the bremsstrahlung continuum is 
especially well sampled, even weak higher temperature com- 
ponents can significantly contribute to the flux. Also, these 
components would lead to a higher average temperature for 
the ICM than if the average cluster temperature were deter- 
mined from the spectrum at softer energies, such as from 0.5- 
lOkeV. 

In the previous section, we found that a two temperature 
model provided a slightly better description than did a sin- 
gle temperature model of joint fits to the XMM-Newton and 
Suzaku data, especially when ignoring the 3% systematic er- 
ror applied to the XMM data. This may indicate that there 
are multiple temperature components in Coma. The multi- 
ple components could occur along the line-of-sight, or in the 
plane of the sky, or locally (the gas might be multiphase). In 
fact, previous temperature maps show that Coma certainly has 
temperature structure w hich is likely associated with mergers 
( Watanabe et al.l 1 1999th Here, we test whether this tempera- 
ture structure alone could reproduce the observed Suzaku PIN 
spectrum of Coma, without any non-thermal emission. From 
the XMM-Newton EPIC-pn mosaic of Coma, we constructed 
a temperature map on a 16 x 16 grid with cell size of 4. '3 on 
a side. Each of the spectra were fit with a single tempera- 
ture APEC model to produce a temperature map that covers 
the Coma mosaic, as shown in Figure [7] We weighted these 
model fits by the PIN spatial sensitivity and combined them. 
This resulting model was compared to the PIN spectrum (Ta- 
ble|2] row labeled "T map "). Note that only the overall normal- 
ization of the T map model was allowed to be fit, to compen- 
sate for a loss of flux due to incomplete coverage of the map 
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FIG. 8. — Suzaku HXD-PIN spectrum (E > 12 keV) and the combined 
XMM spectrum (E < 12 keV) corresponding to the spatial sensitivity of the 
PIN. Shown as solid lines are the combined spectra of the best fit models 
from the temperature map. The thermal model ("Tmap", green) is nearly 
coincident with the data, though falling below it at higher energies. The other 
two components are described in Figuref4] 

across the HXD FOV. Also, the spectral shape and normaliza- 
tion of each of the thermal models was the same as given by 
the XMM-Newton temperature, and each model was simply 
weighted by the average PIN effective area at that position. 

This T map model provides a good fit to the PIN spectrum 
with no adjustable parameters that affect the spectral shape 
(Fig. [8}. While the value of % 2 is worse than for all previous 
fits, it is only slightly worse, which is not unexpected since 
it is not really a fit at all. The fact that XMM-derived tem- 
peratures extrapolated to PIN energies are sufficient to fully 
account for the PIN spectrum further suggests that no non- 
thermal hard X-ray excess has been detected with the PIN, 
especially below 40 keV and probably below 70 keV. Also, 
simulating spectra of similar quality to our data, assuming the 
Tmap model for the underlying source, yields a joint single 
temperature fit consistent with that found from the actual data, 
with T = 8.51 ± 0.06 keV. 

4.3. Systematic Errors in the Spectrum 
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FIG. 9. — Suzaku HXD-PIN spectrum (E > 12 keV) and the combined 
XMM spectrum (E < 12 keV) corresponding to the spatial sensitivity of the 
PIN. The PIN NXB is increased to its 90% confidence limit, which is 2.3% 
for E < 40 keV and 4% for E > 40 keV. The thermal model ("APEC", green) 
is nearly coincident with the data, though falling below it at higher energies. 
The other two components are described in Figure [4] Note that the residuals 
above the model that exist in the previous fits at E > 40 keV have disappeared. 



We explicitly consider the systematic error for 3 quanti- 
ties: the PIN non-X-ray background, the XMM-Suzaku cross- 
normalization factor, and the normalization of the CXB as 
modeled for the PIN spectrum. To test the effect of these sys- 
tematic errors, we vary the relevant quantity up and down by 
our estimate of the 90% systematic error, and evaluate the re- 
sulting change in best-fit model parameters. The detection of 
a non-thermal component cannot be claimed unless it remains 
robust to variations of these quantities within their systematic 
errors. Because the largest error is in the normalization of the 
PIN NXB, we first increase it by 2.3% from 12-40 keV and 
4% from 40-70 keV and repeat the single temperature plus 
non-thermal model fit. The new best-fit IC normalization is 
pushed to zero. Even before considering the effect of other 
systematic errors, from this exercise alone it is clear that we 
do not detect non-thermal emission in the HXD-PIN spec- 
trum, given the current uncertainty in the NXB normaliza- 
tion. This fit, with a temperature for the thermal component 
of 8.33 ±0.06 keV, is shown in Figure [9] In fact, the value 
of x 2 is slightly lower (x 2 = 1672.25 for 1689 dof) than for 
the nominal PIN NXB single temperature fit. Notice that the 
residuals above the model for E > 40 keV, seen in all the pre- 
vious spectral fits, have disappeared. 

Though we cannot claim to detect non-thermal emission, 
we can derive an upper limit to its flux based on joint fits to 
the spectra, including systematic errors in the following way. 
First, for an assumed photon index which we fix, we find the 
nominal normalization N nom and corresponding 90% upper 
bound N^ om of the non-thermal component for a single tem- 
perature plus power law model, allowing the temperature and 
normalization of the thermal component to vary. Then, for 
each systematically uncertain quantity, we set that quantity to 
the limit bounded by the systematic error in the sense that in- 
creases the value of the non-thermal normalization N sys j, and 
we fit for it and its 90% upper bound ,. The statistical 
and systematic errors of the power law normalization are then 
given by 

S^ t =Kl m -N nom , (11) 

and 

^sys,i = ^sys./ — ^noni = ^sys,/ — -Miom i (12) 
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FIG. 10. — Suzaku HXD-PIN spectrum (E > 12 keV) and the combined 
XMM spectrum (E < 12 keV) corresponding to the spatial sensitivity of the 
PIN. All quantities with systematic uncertainties (PIN CXB and NXB, the 
XMM-Suzaku cross calibration) are set to their 90% confidence limit in the 
direction that favors the addition of a power law model component to de- 
scribe the data. The thermal model ("APEC", green) is nearly coincident 
with the data, though falling below it at higher energies. The non-thermal 
model ("Power Law", light blue) is shown for T = 2 at its 90% confidence 
upper limit value. The other two components are described in Figure|4] 

respectively. The final 90% upper limit is then given by 



< = Nnom + ^ S W,> + ■ d3) 

We add each systematic error contribution in quadrature be- 
cause it is unlikely that we chose normalizations for these 3 
quantities such that each one disfavors the detection of non- 
thermal emission in the most severely possible way. 

The upper limits for a range of assumed photon indices is 
provided in Table [3] and in Figure [10] we show, for T = 2, 
the resulting best fit with all 3 systematic quantities set at 
the limit of their 90% confidence range. Over the PIN en- 
ergy band (12-70 keV), the flux is relatively independent of 
photon index. To compare our results to the most recent pre- 
vious detections of non-thermal emission in Coma, we also 
give the upper limit on the non-thermal flux in the 20-80 keV 
band, which is 6.0 x 10~ 12 erg s" 1 cm" 2 for T = 2. This limit 
is inconsist ent with the RXTE (Rephaeli & Gruber 20Q2) and 
BeppoSAX (iFusco-Femiano et alJ 12004) detections by about 
a factor of 2.5, but at the same le vel as the upper limit de- 
rived bv lRossetti & Molendil d2004|) from the BeppoSAX data.. 
If we adopt a 4% systematic error for the PIN NXB instead 
of 2.3% for E < 40 keV, which would better agree with that 
derived from "blank sky" observations, then our upper limit 
increases by 35%. Similarly, if we also increase the CXB and 
XMM-Suzaku cross-normalization to 18% and 10%, respec- 
tively, our upper limit for T = 2 would increase by 50%. In 
either case, our upper limit still excludes the RXTE and Bep- 
poSAX detections, if FOV differences are ignored (see §[5]for 
a more meaningful comparison). 

5. IMPLICATIONS AND DISCUSSION 

After modeling all the known possible contributions to the 
2-70 keV spectrum, simultaneously fitting for the param- 
eters of thermal and non-thermal spectral components, and 
taking into account the systematic uncertainty of the PIN 
NXB, we do not see evidence for IC emission in Coma 
at our level of sensitivity. We therefore derive an upper 
limit to non-thermal, hard X-ray emission through a care- 
ful consideration of the maximum effect of systematic un- 
certainties on our ability to detect a non-thermal signal. 
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TABLE 3 
90% Upper Limits on IC Flux 



Norm. a Flux (12-70 keV) Flux (20-80 keV) BeppoSAX Detection 15 

T (10~ 3 photons keV" 1 cm" 2 s" 1 ) (10~ 12 erg s" 1 cm" 2 ) (10~ 12 erg s~' cnT 2 ) (10~ 12 erg s" 1 cnT 2 ) 



1.0 


0.155 


14.4 


14.9 


1.1 


0.220 


14.2 


14.4 


1.2 


0.311 


14.1 


13.9 


1.3 


0.439 


13.9 


13.4 


1.4 


0.617 


13.8 


12.9 


1.5 


0.860 


13.5 


12.3 


1.6 


1.18 


13.1 


11.6 


1.7 


1.58 


12.4 


10.7 


1.8 


2.04 


11.4 


9.51 


1.9 


2.48 


9.83 


7.98 


2.0 


2.70 


7.64 


6.01 


2.1 


2.70 


5.46 
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11 Normalization of the power-law component for the T+IC model, which is the photon flux at a photon energy 

of 1 keV. ; 

b Flux (20-80 keV), as reported in Fusco-Femiano et al. 1 2004). 



This conservative upper limit is similar to that derived by 
iRossetti & MolendU (2004) a nd is inconsistent wi th claimed 
detectio ns using RXTE dRephaeli & Gruberll2002l) and Bep- 
poSAX dFusco-Femiano et alj|2004l) by approximately a fac- 
tor of 2.5. However, it should be noted that we do not include 
potentially lost emission due to PIN vignetting from any of 
our flux upper limits relative to the larger FOVs of RXTE and 
BeppoSAX, which have collimator FWHM of 1 and 1 .3°, re- 
spectively. 

If IC emission follows the radio synchrotron emission [as 
derived from the point source-subtracted radio image from 
iDeiss et all (119971) 1. as it would for a uniform B field through- 
out the cluster, our upper limits imply a total IC flux 1 .7 - 2 x 
larger would be found inside an RXTE/BeppoSAX-liks FOV. 
We also consider a more detailed spatial distribution for 
the under lying IC emission , deriv ed from the re-acceleration 
model of iBrunetti & Blasl d2005l) . in which relativistic pro- 
tons collide with electrons in the ICM that are then re- 
accelerated by Alfven waves due to cluster mergers. Given 
the radio spectrum of Coma, this model predicts that the 
smaller FOV of the Suzaku HXD-PIN would lead to an un- 
derestimate of the non-thermal flux by a factor of 2-2.5 
(possibly 3 under extreme circumstances). If this model for 
the spatial distribution of the non-thermal emission is cor- 
rect, the n our upper limit is just co nsist ent with the measure- 
ments of Repha eli & Grub er (2002) and Fusco-Femian o et all 
(2004). (However, the IBrunetti & Blasil d2005h model actu- 
ally predicts a non-thermal flux considerably below the Bep- 
poSAX and RXTE detections.) Because any spatial variation 
of the magnetic field strength is unknown, a direct compar- 
ison between these missions is not possible with any preci- 
sion. Under the reasonable assumption that B decreases with 
radius, our upper limit will be 2 x larger, so our result can- 
not definitively rule out the detections discussed above. Re- 
gardless of this issue, the upper limit is approximately the 
sa me as or slightly higher th an the upper limit range found 
by Rossetti & Molendi (2004). 

However, the BeppoSAX PDS measur es a 20-80 keV flu x 
for the Crab of 1.23 x 10" 8 erg s" 1 cm" 2 dKirsch et al.ll2005l) . 
while the Suzaku PIN flux over this energy range is 1 .56 x 
10~ 8 erg s -1 cm -2 , after applying the 13.2% correction to bring 



the PIN spectrum into agreement with the XIS fluxes (Ishida 
et al., Suzaku Memo 2007-1 1 14 ). This 21% flux difference 
implies our upper limit would be 4.7 x 10~ 12 erg s" 1 cm" 2 on 
the BeppoSA X scale, which is on t he low er end of the range 
estimated by Rossetti & Molendi (2004J). Also, even if we 
only detect one-third of the total emission observed by the 
BeppoSAX PDS, we would just barely exclude the nominal 
value of the lFusco-Femiano et al] d2004l) measurement. 

Assuming the differing measurements of non-thermal emis- 
sion are not due to the IC radiation having a larger extent, 
what might be the cause of this discrepency? While it could 
be explained by a greater point source contamination at hard 
energies for the RXTE and BeppoSAX missions due to their 
larger FOVs, most likely we differ in our results because of 
different considerations of the thermal gas. For both detec- 
tions, the gas temperature was found to be lower than our 
nominal value of 8.45 keV. Fixing the gas temperature to their 
assumed values in our fits yields a T = 2.0 non-thermal com- 
ponent significance > 4er, without including systematic ef- 
fects, for T = 7.67 keV (RXTE) and T = 8.2 keV (BeppoSAX); 
however, these fits are poor relative to fits in which the tem- 
perature is a free parameter. Though these temperatures differ 
from our best-fit value by only a few percent, the exponential 
decline of bremsstrahlung continuum at high energies ampli- 
fies even small differences. The lower measurements of the 
ICM temperature appear not to be due to the inclusion of data 
at low energies (E < 1 keV), which can bias average tempera- 
ture estimates low. Most likely, the larger FOVs of RXTE and 
BeppoSAX allowed the inclusion of emission from more cool 
gas in the cluster outskirts than was observed by Suzaku. This 
emission would serve to lower the average observed temper- 
ature, which is primarily determined from emission at lower 
energies (E < 10 keV). But, as evidenced by the temperature 
map in Figure|7] a distribution of higher-than-average temper- 
ature regions can effectively increase the avera ge gas tempera- 
ture ob served at high energies, as first seen by iNeumann et alJ 
(2003). 

We take the reasonably good agreement between the ther- 
mal models derived from the XMM temperature map and the 

14 http://www.astro.isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo- 
2007-1 l.pdf 
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PIN spectrum to mean that we essentially only detect thermal 
emission from Coma out to 70 keV. This result is fully consis- 
tent with re cent INTEGRAL det ections of extended hard X-ray 
emission. Ren aud et al.l (12006b performed a point-by-point 
spectral comparison between XMM-derived and INTEGRAL- 
derived temperatures and found that they fol l owed a strict 
one-to-one correlation. Similarly, Eck ert et ail (|2007a) char- 
acterized a surface brightness excess relative to the XMM data, 
which they found to be best described by extended hot, ther- 
mal emission ataT~12±2 keV. This excess coincides with 
the hotter temperatures (T ~ 10- 11 keV) to the west of the 
PIN pointing center in Figure|7] 

From our upper limit on the flux of IC emission, we can 
derive a lo wer limit on the average magn etic field strength B 
as shown by Harris & Romanishin ( 1974). Equation (O refers 
to the total energy emitted for one electron; it is more useful to 
consider the ratio of monochromatic fluxes Fr(vr), Fx(vx) at 
frequencies vr, vx, for a power law distribution of electrons, 
from which we can derive an expression for the magnetic field 
as 

B = C(p)(l + z)^ +1) [^-j rA , (14) 

where p is the index of the electron distribution N(E) oc 
E~ p and is related to the spectral index a (F„ oc v~ a ) by 
p = 2a + 1. The value of the proportionality consta nt C(p) 
can be found from the ratio of th e synchrotron flux (lLongai: 
[1991 eqn. 18.49) to the IC flux dRybicki & Lightmanlll97 " 
eqn. 7.31). Assuming that the electron energy distribution 
does not turn over significantly at low energies and that a = 1, 
we find B > 0.15 uG. This limit is still be low the equipartition 
value of 0.5 fiG dGiovannini et al.l [T993). and it is well below 
the line-of-sight estimates of a few ij,G derived fro m Faraday 
rotation measure (RM) studies dFeretti et alJl!995 l). Note that 
the Faraday RM magnetic field estimates are sensitive to the 
B field geometry and may imply a field strength larger than 
the volume-averaged value if B is pref erentially aligne d along 
filamentary structures on small scales (Petrosian 2001). Also, 
we are unable to put interesting constraints on the relativis- 
tic energy budget of Coma, since our lower limit includes the 
equipartition estimate of B, which defines the minimum en- 
ergy in relativistic components and would not imply a signif- 
icant contribution to the energy budget of Coma. 

It has been not ed that the hard excess detected by 
lEckert et alj d2007al) also corr esponds to the peak in the point 
source-subtracted image from lDeiss et al.l (H997), potentially 
indicating that the hard emission could in fact be non-thermal 



in origin (Ecke rtet alj|2007bi) . In fact, we suspect that this 
peak, which appears tantalizingly like a small radio relic, is 
not a true feature of the halo, but instead is the result of im- 
perfect source subtraction. Due to the large beam size used 
to create the diffuse radio image, extended emission from ra- 
dio galaxies might not have been properly subtracted using a 
point source list. We point out that the strongest radio source 
(1256+282 or 5C 4.81, centered on NGC 4869) in Coma is 
near this position, is a head-tail rad io galaxy with a steepe r 
spectral index than rest of the halo dGiovannini et alJll993h. 
and that its tail extends to the west (O'D ea & Owenlll985L 
see Fig. l(j)) and turns north ( Vent uri et al.| [l990, see Fig. 2). 
Subtracting a point source from this morphology would leave 
a residual very similar to that in the Deiss image. Therefore, 
any relation between the location of hard emission and this 
radio feature probably should be regarded as coincidence. 

As the calibration of the NXB model improves, constraints 
on IC emission in Coma will tighten, possibly leading to a 
detection. The uncertainty in the current background model, 
"bgd_d," is more than a factor of 2 lower than the original 
model. However, the existence of non-thermal emission in 
the Coma cluster may have to be determined by future mis- 
sions — particularly, those missions with hard X-ray imag- 
ing capabilities like NuStar 15 , Astro-H (previously NeXT) 16 , 
and Simbol-X 17 . If the IC emission is localized, then our 
joint fitting-technique can be used for many much smaller re- 
gions where temperature-mixing will be less significant and 
the IC component will be relatively stronger. Also, the B field 
strength can be derived spatially across a cluster, yielding a 
better estimate of the possible pressure support provided by 
relativistic components in the ICM, which could modify mass 
estimates that depend on the hydrostatic equilibrium state of 
the cluster gas. 
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